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Refined infra-red magneto-transmission experiments have been performed in magnetic fields B up 
to 35 T on a series of multi-layer epitaxial graphene samples. Following the main optical transition 
involving the n = 0 Landau level, we observe a new absorption transition increasing in intensity 
with magnetic fields B > 26T. Our analysis shows that this is a signature of the breaking of the 
SU(4) symmetry of the n = 0 LL. Using a quantitative model, we show that the only symmetry 
breaking scheme consistent with our experiments is a charge density wave (CDW). 


I. INTRODUCTION 

In multicomponent quantum Hall systems, interaction 
effects lead to a rich variety of broken symmetry ground 
states. In graphene, the spin and valley degrees of free¬ 
dom of the lowest Landau level (LL) form an SU(4) 
symmetric quartet. Refined transport experiments have 
shown evidence of a broken symmetry state , but 
there is no clear consensus on its nature mi , and how 
it is affected by different substrates and disorder. While 
the spin deCTee of freedom has been probed in tilted mag¬ 
netic fields [l|,0)S0, we show here that the valley degree 
of freedom can be accessed by examining the signatures 
of optical phonons in magneto-transmission spectra. In 
this paper, we show that our observation of a new absorp¬ 
tion transition supports the existence of a charge density 
wave (CDW) in our epitaxial graphene samples. 

In a quasiparticle picture, charge carriers in graphene 
are characterized by a Dirac-like spectrum around the 
K and K' equivalent points (“valley”) of the Brillouin 
zone of the hexagonal crystal lattice. As a consequence, 
the application of a magnetic field B perpendicular to 
the plane of the structure splits the electronic levels into 
Landau levels (LL) indexed by n, with specific energies 
En = sgn{n)vF^y2ehB\n\ where n are integers includ¬ 
ing 0 (r^F being the Fermi velocity). In this paper, we 
are concerned with how a broken symmetry phase can 
be observed in infra-red magneto-optical transitions in¬ 
volving the n = 0 LL (i.e., transitions from n = — 1 to 
n = 0 or from n = 0 to n = 1 equivalent to a cyclotron 
resonance (GR) transition in the quantum limit) with an 
energy Eqi = vpV^eKB [l^. Our previous work 0 has 
reported on the magnetic field dependence of this tran¬ 
sition revealing its interaction with the AT-phonon. Be¬ 
sides this specific interaction, we observed that the basic 
broadening 701 (H) oc y/B of the transition had an ad¬ 


ditional component proportional to B in contrast to all 
theoretical models [l^. This could be already a sign of 
the breaking of the valley degeneracy. 

In the present work, we use the F-phonon at the Bril¬ 
louin zone center as a probe of the valley symmetry 
breaking. In the absence of valley symmetry breaking, 
the F-phonon does not affect the infra-red absorption 
spectrum because the electron-phonon matrix elements 
are of opposite signs for the K and K' valleys [l^ . How¬ 
ever, one expects to see signs of valley symmetry break¬ 
ing when the energy Eqi{B) is larger than that of the 
optical F-phonon {tujjp = 0.196eV). It turns out, indeed, 
that when that condition is reached, a new optical tran¬ 
sition develops at an energy higher than the main line 
(Fig. [T]). We interpret this as a signature of the breaking 
of the SU(4) symmetry. A model has been established 
to reproduce these findings and applied to the different 
phases which have been proposed. 

In Sec. im we discuss these experimental observations 
and methods in more detail. We first interpret our ex¬ 
perimental findings within a simplified model for valley 
symmetry breaking in Sec. IHII before deriving a more 
complete Hamiltonian in Sec. llVI and calculating the op¬ 
tical conductivity in various broken symmetry phases in 
Sec. 13 Finally, a comparison between experiment and 
theory is presented in Sec. IVII followed by conclusions in 

Sec, iwn 


II. EXPERIMENTAL OBSERVATIONS AND 
METHODS 

In our experiment, precise infra-red transmission 
measurements were performed on multi-layer epitaxial 
graphene samples, at 1.8 K, under magnetic fields up 
to 35 T. The light (provided and analyzed by a Fourier 
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FIG. 1: (color online). Left panel: Evolution of the Eoi tran¬ 
sition for different values of the magnetic field beyond 27 T. 
Experimental transmission data for sample S5 (open circles) 
is compared with calculated transmission spectra (red lines), 
for different values of magnetic field, using the proposed CDW 
model. Right panel: 1 - (relative transmission) measured in 
the experiment and calculated (red lines) for the CDW phase, 
for magnetic fields B = 27, 29, 31, 33, 35 T. Deconvolution 
of the experimental spectra into two Lorentzians is shown in 
blue dashed and dotted lines. 


transform spectrometer) was delivered to the sample by 
means of light-pipe optics. All experiments were per¬ 
formed with nonpolarized light, in the Faraday geom¬ 
etry with the wave vector of the incoming light paral¬ 
lel to the magnetic field direction and perpendicular to 
the plane of the samples. A Si bolometer was placed 
directly beneath the sample to detect the transmitted 
radiation. The response of this bolometer is strongly de¬ 
pendent on the magnetic field. Therefore, in order to 
measure the absolute transmission TA {B,uj), we used a 
sample-rotating holder and measure for each value of B 
a reference spectrum through a hole. These spectra are 
normalized in turn with respect to TA(0, cu) to obtain 
a relative transmission spectrum TR(i3,a;) which only 
displays the magnetic field dependent features. Those 
spectra are presented in Fig. [2l 

The samples were grown 17[ on the C-terminated sur¬ 
face of SiC and display the characteristic transmission 
spectra of isolated graphene monolayers that arise from 
rotational stacking of the sheets [l^. The thickness d of 
the SiC substrate has to be reduced significantly in order 
to minimize the very strong double-phonon absorption 
of SiC in the energy range of interest. In the first se¬ 
ries d was reduced to 60^m and related samples have 
been used to perform the experiments reported earlier 
[ 3 . One of them, named S4, was used to compare the 
data with those obtained on sample S5 from a new se¬ 
ries where the thickness d was further reduced down to 
32/iTO. We compare in Fig.[2]the transmission spectra, at 
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FIG. 2: (color online). Top: Relative transmission spectra 
of sample S5, for different magnetic field values up to 35 T. 
Bottom: Evolution of the Eqi transition for different values 
of the magnetic field between 27 T and 34 T, for samples S4 
(open dots) and S5 (full lines). 


high fields, for samples S4 and S5. Technically speaking, 
the optical response of both samples is almost the same, 
showing that they have a similar number of active layers. 

Taking into account all layer dielectric properties of 
each sample in a multi-layer dielectric model, we de¬ 
termine the effective number, N^ff, of graphene sheets 
with their respective carrier density (see Appendix). For 
samples S4 and S5, we have found that N^ff = 7 with 
carrier densities {5.5, 3.2,1.8,1.1, 0.5, 0.1,0.1}10^^cm“^ 
and {6.0, 3.5, 2.0,1.2, 0.5, 0.1, 0.1}10^^cm“^ respectively. 
These carrier densities are fixed for each sample. 

The transmission spectra of sample S5 at high mag¬ 
netic fields are displayed in Fig. [TJ We observe a new 














3 


transition occurring at an energy higher than that of the 
CR line, growing in intensity when increasing the mag¬ 
netic held. This behavior cannot be explained without 
breaking the SU(4) symmetry in graphene. In order to 
characterize more clearly these hndings, one can treat the 
data, as a Hrst step and in a very rough way, extracting 
from the transmission data the real part of the effective 
diagonal component of the conductivity axxi^^B) [ 2 nj |. 
We have deconvoluted this result with two Lorentzians of 
equal width, extracting the evolution of the two extrema 
with the magnetic held. The resulting energies are dis¬ 
played in Fig. [3] (top panel) for samples S4 and S5. 

Though the procedure adopted at this initial level is 
quite rough, it provides important information: (i) The 
evolution of the lower energy line varies at low helds like 
(function F2{B) in Fig. ED with a coefficient pro¬ 
portional to the Fermi velocity vp and ends at higher 
helds with a similar dependence (function F1(R)) but 
with a smaller value of vp which is, by itself, a sign of 
some interaction occurring at an energy close to that of 
the F-phonon; (ii) The second component of the decon¬ 
volution always appears at energies larger than that of 
the F-phonon; (hi) In principle, in the SU(4) symmet¬ 
ric picture, it is not possible to explain the occurrence 
of an additional transition, growing in intensity with B, 
at higher energies than the main transition line; (iv) It 
is therefore clear that the F-phonon plays a crucial role 
though it should not, indicative that the SU(4) symmetry 
is broken. Using these observations, we now have some 
guidelines to develop a theory which can explain quanti¬ 
tatively the experimental observations. In addition, we 
note that results for samples S4 and S5 are quite similar 
within the experimental errors. Knowing that the ac¬ 
tive layers which contribute to the Eqi transition should 
have a filling factor v < 2 ( zz = A^s<I>o/i?, $0 being the 
flux quantum and Ng the carrier density) and that, in 
samples S4 and S5, the carrier density for active layers 
do not have the same sequence, the physical mechanisms 
describing the experimental findings should not be very 
dependent on the doping of active layers. This is indeed 
the case as discussed below. 


III. SIMPLIFIED MODEL FOR VALLEY 
SYMMETRY BREAKING 

To illustrate how the electron-phonon interaction and 
the valley symmetry breaking give rise to the observed 
features in the transmission spectrum, we first introduce 
a simplified model for the interaction of the F phonon 
with the Uoi excitation, before discussing the full SU(4) 
calculation. The simplified model provides a minimal de¬ 
scription of the valley symmetry breaking by neglecting 
the spin degree of freedom in the n = 0 LL. We assume 
that K and K' sublevels of the n = 0 LL are separated in 
energy by Ay, and have different filling factors pk and 




FIG. 3: (color online). Top: CR energies, resulting from the 
deconvolution of the experimental transmission traces, as a 
function of y/B for samples S4 (full dots) and S5 (open dots). 
Error bars are similar for both samples. The full lines are 
a linear fit of the data for the low energy transition at low 
(F2(i3)) and high (Fl(i3)) field. Bottom: Comparison of the 
variation of experimental CR energies for sample S5 (open 
dots), with that calculated for the CDW phase (open squares), 
as a function of y/B. The size of the open squares mimics the 
relative oscillator strength of the optical transition. 


pk'- Considering just the n = 0 to n = 1 transitions, the 
interaction with the F phonon is captured by the Hamil¬ 
tonian (in the basis of creating an electronic excitation in 
K, electronic excitation in AT', and a F phonon, in that 
order) 


^Eqi — Ay/2 0 gphy/^ \ 

H~\ 0 Uoi-l-Ay/2 —gphy/PK' 1 ■ (1) 

gphy/Ep —gphy/PK' ) 


where gph characterizes the electron-F-phonon interac¬ 
tion. The optical conductivity is calculated using the 
Green’s function formalism introduced by Toyozawa 211. 
The diagonal component of the conductivity is: 
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axx{f^) = —lmM^G{fujj)Mx ( 2 ) 

(jJ 

where the Green’s function is G = {huj — H — iri)~^, with 
77 —>• 0+ (see Sec. IIVBT) . The optical matrix elements for 
the simplified model are \/vk'i O) . The sim¬ 

plified model explains the splitting of the main transition 
line in the two limits Eqi rs Hujr and Eqi ^ hur- In the 
absence of valley splitting (Ay = 0 and vk = vr'), the 
eigenstates of the pure electronic part of H in Eq. [T] are 
valley-symmetric and valley-antisymmetric combinations 

of Eqi transitions, i.e. ^ rCo,k + A.k-c^,k) |GS) 

and \ rCq^k — cj riCq^ri'^ |GS), respectively, where 

IGS) denotes the ground state and ^ are creation op¬ 
erators at LL n and valley K. The valley-symmetric com¬ 
bination is infra-red active but does not interact with 
the r phonon, while the valley-antisymmetric combina¬ 
tion is infra-red inactive and interacts with the T phonon. 
The symmetry breaking valley splitting term Ay allows 
both eigenmodes to interact with the E phonon while 
remaining infra-red active, inducing a splitting of the 
main transmission line in the vicinity of the T phonon 
frequency. Away from the T phonon frequency {Eqi ^ 
hujr), Eqi transitions at K and K' interact weakly with 
the phonon and the splitting of the main transmission 
line is controlled directly by the energy difference Ay. 


IV. THEORY OF MAGNETO-PHONON 
RESONANCE IN THE PRESENCE OF SU(4) 
SYMMETRY BREAKING 


We reintroduce the spin degree of freedom and the 
n = —lton = 0 transitions in order to obtain a quantita¬ 
tive understanding of the experiment. We consider differ¬ 
ent theoretical models of the n = 0 LL SU(4) symmetry 
breaking, taking into account the effects of y 0 and dis¬ 
order by introducing Gaussian broadening into a mean 
field theory (Sec. IIVAD . Different symmetry-breaking 
phases are represented in the mean field theory by dif¬ 
ferent orderings and filling factors of the four sublevels 
of the n = 0 LL. We consider four candidate symmetry¬ 
breaking phases that have been proposed in the literature 
0: Ferromagnetic(F), Charge Density Wave (CDW), 
Canted Antiferromagnetic (CAF) and Kekule-distortion 
(KD), and calculate the optical conductivity using Eq. [2] 
with the appropriate Hamiltonian H for each phase. 
Treating these phases on the same footing (detailed in 
Sec. El, we find that each phase results in characteristic 
features in the evolution of the transmission spectrum as 
a function of the magnetic field. By examining the inten¬ 
sities and positions of the transmission lines, we identify 
the symmetry broken phase in the samples used in our 
experiment as the CDW type 271. l28|. 


A. Description of the ground state 

We assume that the ground state is a single Slater de¬ 
terminant of the form: 


4 Nj 

1^^)=n n(3) 

j=l mj=l 

the index j runs over the 4-dimensional spin/valley space 
and VTij describe the ’’guiding center” degree of free¬ 
dom. The state is represented by the wavefunc- 

tion where is a four-component spinor and 

is the orbital part of the wavefunction. These 
wavefunctions belong to the n = 0 landau level (LL) of 
graphene. The occupation numbers Nj count the number 
of j states that are occupied in this ground state. 

There are different models proposed to describe the 
symmetry-broken phase of graphene which have been re¬ 
viewed by Kharitonov . For a given model, we assume 
that the system is polarized along a certain direction in 
/-space. For instance, with increasing order of energies, 
j = 1,2,3,4 corresponds to (AT' l,K 1) in 

the charge density wave (CDW) phase. The remaining 
degrees of freedom, (j)mji^) and Nj, are treated as vari¬ 
ational parameters, subject to the constraint 7Vi -|- V 2 -I- 
A 3 + A ^4 = N. We minimize the energy of the ground 
state Eqq = {GS\E{) -t- -t- E(iisQxder\GS'). Here Hq 

is the single part of the Hamiltonian without disorder, 
iLe_e the interaction term and Eldisorder the disorder po¬ 
tential. Because we assume a single Slater determinant, 
we can apply mean-field theory and obtain single-particle 
energy levels Ej^rnj (The origin of the energies is taken to 
be at the energy of the n = 0 LL of the non-interacting 
system). 

In a system with finite disorder, the energy levels Ej^rnj 
are clustered about mean values Ej = avg Ej^ruj ■ We re- 

move the rrij degrees of freedom by replacing the energy 
levels Ej^ruj by broadened energy levels centered at Ej. 
There is a Fermi level Er which fixes the occupation 
numbers Nj when the graphene layer is doped with a 
total filling factor v. Assuming the broadening to be of 
Gaussian type with a width 70 the Fermi level is deter¬ 
mined by solving the following equation: 


V = 


y;Eif( 


Er — Ej 


(4) 


from which one can calculate the individual filling factors 
Vj = (l-|-Erf(-^^^^))/2 for each level Ej. These Ej will 
be used, later on, as fitting parameters dependent on the 
broken-symmetry phase under consideration. Note that 
in this approach all optical transitions to or from the 
n = 0 LL are allowed. 
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B. Description of the optical transitions 

We first consider the transitions from the n = 0 LL 
to n = 1 LL. The Hamiltonian of the magneto-excitons, 


including their interaction with the T-phonon, denoted 
iLo (reminding that it describes the optical transitions 
allowed in the polarization), is: 


(hjjjQi — El 0 

0 fvjJm — E 2 


Hn = 


V 9 iV^ 


0 

0 

92 


0 0 5iv^\ 

0 0 52v^ 

HlOq! — E 3 0 <?3-\Ah 

0 huoi - E 4 , 

93V^ 9X\/Ei f^ph) 


V 


(5) 


where fia;oi is the energy of the Eqi transition from the 
n = 0 LL to n = 1 LL in the absence of interactions 
and hiOph that of the T-phonon. This Hamiltonian de¬ 
scribes the excitations from the 4 sublevels {Ej^j = 1..4} 
of the n = 0 LL to the n = 1 LL. The matrix ele¬ 
ments j = 1..4} respectively describe their interac¬ 
tion with the T-phonon, and is dependent on the wave- 


function character of the 4 sublevels (i.e., dependent 
on the broken-symmetry phase). In general, gj oc y/B 
[H, with a prefactor dependent on j and the broken- 
symmetry phase. 

Similarly, the Hamiltonian describing the magneto- 
excitons for the transitions from the n = —1 LL to the 
n = 0 LL (allowed in the a~ polarization) is written as: 


Ho = 


/ huJoi + El 0 

0 Eoqi E 2 
0 0 

0 0 

\5lVl - S 2 VI - V2 


0 0 

0 0 

hiUQl + 0 

0 Hujqi E/^ 
ffaVl - V 3 gl^/l - 


gi Vl - i/A 
g2\ /l ~ V 2 
g3\ /l - V 3 
giy/l - VA 

J 


The total Hamiltonian H describing the magneto- 
excitons is therefore: 


H = 


(Ho 

I 0 



(7) 


CTxx(w) = -mJ.g.Mj, 

OJ 

a^yiuj) = -m^.g.m; 

UJ ^ 


( 6 ) 


( 8 ) 


We will also need to introduce the optical matrix ele¬ 
ments Mx and My for corresponding transitions. These 
matrix elements depend on the ground state under con¬ 
sideration. In CDW case they are (see Sec. IV Ap : 
Mx /vq = {yd ^, \fv 2 , \fvi-, y/vi, 0, -y/1 - Vi, -y/1 - V 2 , 
— y/1 — 1 ^ 3 , — ^/l ~ ^ 4 j 0} and My = iMx- For a different 
scenario, the optical matrix elements will be transformed 
to a different basis, as will be detailed in Sec. El 

The Green’s function for the magneto-excitons, is ob¬ 
tained as G = {{hu! + i 7 oi).I — H)~^ (where I is the unit 
matrix and 701 the broadening of the Eqi transition). 
This allows us to calculate the different components of 
the conductivity: 


V. OPTICAL CONDUCTIVITY IN THE 
DIFFERENT PHASES 

Here, we calculate the optical conductivity for the dif¬ 
ferent symmetry broken phases, using Eq. [8l 


A. Charge density wave (CDW) phase 

The CDW phase is characterized, at filling factor 
= 0, by two electronic LL full in one valley (say K' 
for instance) and two LL empty in the other valley (AT). 
We will introduce a valley asymmetry ZW mainly deter¬ 
mined by electron-electron interactions (9| and a Zeeman 
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splitting As- Therefore the sequence of sublevels take In this case, the parameters governing the electron- 
the following form: Tphonon interaction gi, g 2 on one hand and ( 73 , g^ on 

the other hand are of opposite sign. That is, 

Ei{K' t) = -Av/2-As/2 
i:;2(i^';) = -Av/2 + A5/2 
E 3 {K\) = Av/2-As/2 

- E4(K ;) - Av/2 + As/2 - 1 


{GScdw + rphonon|i7e-p/t'I'J^g^^Tif5^o|G'S'c_Dw) = gphlV2 
{GScdw + rphonon|i/e_p?i\I']^,^^^v[/^,j,_Q|G'S'c_Dw) = —gphlV2 
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FIG. 4: (color online). CDW phase with Av oc B-. evolution 
of the Oxx component of the conductivity for different values 
of the magnetic field between 25 T and 35 T for a) a carrier 
density Ns = 0.5 x 10^^cm“^ and b) = 1 x 10^^cm“^. In 
both cases gph = 2.3 x y/ B[T] meV. 


where ^ „ is the creation operator for electrons in val¬ 
ley A", spin s, Landau level n. On the other hand, the 
electron-light interaction (which determines M) has the 
same sign at both valleys. 


{GScdw + piioton\He-Hght'^Ks^i'^Ks,o\GScDw) = 1 

{GScdw + photon|iLe_i ight^^K's K'8g\GScDw) = 1 

( 11 ) 

The results obtained for this phase are presented in 
Fig.m assuming Ay proportional to B, for two extreme 
values of the carrier density. The electron-phonon cou¬ 
pling was taken to be gph = 2.3 x B[T] meV, which 
agrees with density functional theory (DFT) calculations 
and experiments The results are not very 

dependent on Ng. The value of As = 0.15 meV B corre¬ 
sponds to a g-factor of 2.6 to be compared with 2.7 ± 0.2 
reported in 22| . The splitting of the transition is directly 
governed by the amplitude of Ay whereas the introduc¬ 
tion of As modifies only the relative amplitude of the 
two transitions. In all cases both Ay and gph need to be 
finite to observe the effect. We finally note that, in this 
case. Ay > As in coherence with the assumption made 


in our previous work 


14l. 
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FIG. 5: (color online). GDW phase with Ay oc y/B: evolution 
of the component of the conductivity for different values 
of the magnetic field between 25 T and 35 T for a) a carrier 
density Ns = 0.5 x 10^^cm“^ and b) As = 1 x 10^^cm~^. In 
both cases gph = 2.3 x ^B[T] meV. 

However there is no clear consensus about the field de¬ 
pendence on Ay [^. Therefore one can alternatively as¬ 
sume that Ay is proportional to y/B. The corresponding 
results are displayed in Fig. keeping all other param¬ 
eters fixed. We obtained essentially the same results as 
in Fig. [4l Within the experimental errors we will not 
be able to differentiate between the two magnetic field 
variations of Ay. 

The CDW state is compatible with the experimental 
results as we will see below. 


B. Kekule-distortion (KD) phase 

In this phase [^, the K and K' valleys hybridize into 

linear combinations AT, AT'. At zz = 0, both spin f and 

spin ], electrons occupy one of these valley-combinations, 

say AT. The y = 0 ground state for the KD phase is 

fk— 10). Therefore, the ’’natural” basis for this 
ift.o iC4,o' ' ’ 

phase, where the density matrix is diagonal, is {K X,K I 
,K' 1} in contrast to the basis {K f, AT ],,K' f 

, K' 4,} used in the CDW phase. Therefore the sequence 
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of sublevels take the following form: 


j = 1 : t 

J = 2 -.Ki 

J = 3 : IT t 
j = A: Wi 


( 12 ) 


The transformation rules for the operators in this basis 
are 


if,. 


K' ,s,n 




(13) 


where 4'— is the creation operator for electrons in val- 

K,s,n 

ley state K, spin s, Landau level n. Making use of this 
change of basis (Eq. [T31) and Eq. [TU] and Eq. [TTl we de¬ 
rive that the electron-light matrix elements do not change 
with respect to the CDW phase, and the electron-phonon 
matrix elements g vanish by symmetry. That is, 


A/g = 5.0 10” cm'^ 



Energy (eV) 


FIG. 6: (color online). Ferromagnetic phase with Ay oc B: 
evolution of the component of the conductivity for differ¬ 
ent values of the magnetic held between 25 T and 35 T for a 
carrier density Ns = 5.0 x 10^^cm“^ and Qph = 2.3 meV/T*^^^. 


{GSkd + Tphonon\He-ph'i’^^ ^Ai^^^qIGSkd) = 0 

{GSkD + photon|i7e_iig;itr|/|^^ ^Al-i^^ g\GSKD) = 1 

_ ’ ’ (14) 

and the same for K'. The structure of the Hamilto¬ 
nian (Eql?!) becomes only diagonal and no splitting is 
observed when calculating the conductivity. Therefore 
the KD phase does not explain the experimental results. 


Ay oc B. The conductivity does not show any significant 
splitting of the main line. In fact there is an eigenvalue 
of the corresponding Hamiltonian larger than that of the 
main line but it remains optically inactive. Therefore 
here also, the E phase does not explain the experimental 
results. 

D. Canted anti-ferromagnetic (CAF) phase 


C. Ferromagnetic (F) phase 

In the F phase the ground state, at filling factor 

= 0, is composed in both valleys K and K' of a single 
full LL with the same spin. In analogy with the CDW 
phase, we will introduce a valley asymmetry Ay and a 
Zeeman splitting Ag. Therefore the sequence of energy 
levels take the following form: 

El (A' ;) = -Ay/2-As/2 
E 2 (E;) = Ay/2-As/2 
E 3 (E't) =-Ay/2 + As/2 ^ ’ 

Ei{K t) = Ay/2 + As/2 

Note that, in this case, Ag should be larger than Ay 
to preserve the ferromagnetic nature of the state. In 
the present case the parameters governing the electron- 
Tphonon interaction (Ea lbIBI) gi, gs on one hand and 52 , 
gi in the other hand are of opposite sign. 

The results are displayed in Fig. [S] where we have taken 
for Ag the same evolution that in the CDW phase and 


The CAF phase for the ground state is described by a 
spin in direction 9k in valley K and a spin in direction 
6 k' in valley K'. (The directions 9k and 9k' are in 
general not opposite to each other except in the special 
case of the anti-ferromagnetic phase). The direction 9k 
is oriented at an angle 9 relative to the magnetic field 
B and the direction 9k' at an angle —9 with respect to 
it. (In the anti-ferromagnetic phase, 9 = 7 r/ 2 ). Here the 
Zeeman splitting should vary like Ag oc cos 9 and if 9 is 
close to 7r/2 this term should not play a dominant role. 
We choose the following order of states: 

Ei(E,0iy) = -Ai/2-A 2/2 
E2(E',0iyO = -Ai/2 + A2/2 
E3(E,7r + 0iy) = Ai/2-A 2/2 
E4(A:',7r-b 01^0 = Ai/2 + A 2/2 

where the introduction of Ai reflects the CAF pattern of 
spin. We assume in addition that the asymmetry between 
valleys is reflected by A 2 (favoring here the K valley). 
To preserve the CAF phase A 2 should be smaller than 




Ai. Similar to the F phase, the parameters governing 
the electron-Fphonon interaction (Ea l5l6l) gi, on one 
hand and 52 , in the other hand are of opposite sign. 
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FIG. 7: (color online). CAF phase with Ai oc %/B and A 2 oc 
'/B: evolution of the component of the conductivity for 
different values of the magnetic field between 25 T and 35 T 
for a) a carrier density Nb = 0.5 x 10^^cm“^ and b) Ns = 
1 X 10^^cm“^. In both cases Qph = 2.3 meV/T^^^. 

The results are displayed in Fig.[3where we have taken, 
Ai and A 2 proportional to a/S- The results are not 
very dependent on the carrier concentration. We observe 
indeed a splitting of the transition when both A 2 and gph 
are different from zero : in fact the splitting is governed 
by A 2 . In the present case we do not have, a priori, a 
guide for choosing the values of Ai and A 2 . In order to 
be consistent with experimental results, we have taken 
for Ai a value which provides an upper transition energy 
close to that observed. 

However the evolution of the spectra does not reflect 
the experimental observations: whatever is the choice of 
parameters, the intensity of the high energy transition 
never reaches that of the main transition in contrast to 
the CDW phase where it should become dominant at 
fields higher than 35 T. This is discussed further in the 
next section. 


VI. COMPARISON OF EXPERIMENT AND 
THEORY 

The KD, F and CAF phases result in transmission 
spectra incompatible with experiment (Fig. [ 8 ]). In the 
KD phase, electrons occupy linear combinations of the 
K and K' valleys; the electron-phonon matrix elements 
vanish by symmetry, resulting in a single transmission 
line. For the F phase, the occupancy of the K and K' 
valleys are almost equal (Sec. IV Cp . and there is no sig¬ 
nificant splitting of the main transmission line (Fig. [5]) 
. Similarly, the calculated transmission spectra for the 
CAF phase show a second CR line of much lower inten¬ 
sity than the main CR line. Deconvoluting these spectra 
with two Lorentzians we hnd a ratio of the CR weights 
of 0.9 ± 0.05 for the experiment, to be compared to the 
value 0.9 for the CDW phase and 0.4 for the CAF phase. 
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FIG. 8: (color online). Real part of the conductivity measured 
in the experiment (open dots) and calculated with the CDW, 
CAF, F and KD phases (red continuous lines), for magnetic 
fields B = 31T (left panel) and 35r (right panel). The De- 
convolution of the model into two Lorentzians is shown for 
the CDW and CAF phases (blue dashed and dotted lines). 


Despite the introduction of valley asymmetry into the 
CAF phase, we find that it cannot explain the observed 
evolution of CR energies in the experiment. 

The CDW phase has unequal occupation numbers of 
the n = 0 LL at the K and K' valleys, corresponding to 
a density modulation of the graphene A and B sublat¬ 
tices in real space. Unlike the ideal disorder-free CDW 
discussed in [9|, both K and K' valleys have non-zero 
occupation number in our calculation, due to disorder- 
induced broadening. Nevertheless, the mechanism giving 
rise to the splitting of the Eqi transmission line remains 
essentially the same as illustrated by the simple model 
(Eq. [T]) above. 

The spin and valley splittings Ag and Ay parameterize 
our model for the CDW phase and determine the Filing 
factors that enter the Hamiltonian and the optical matrix 
elements. We fix Ag usiM the experimental graphene g- 
factor measured in Ref. [22| . We treat Ay as a fitting pa¬ 
rameter, obtaining Ay = 0.24 x H[T]meV. In our calcu¬ 
lations, we use gph = 2.3 x \/i3[T]meV, which is in good 
agreement with density functional theory (DFT) calcula¬ 
tions [ 2 ^ and experiments [H-ii]. We take the position 
of the Eqi transition line and its broadening to be given 
by VF = I-Ol X 10®ms“^ and 701 [meV] =3-1- 0.8\/i3[T] 
respectively, consistent with their measured values at low 
magnetic fields away from the F-phonon frequency. For 
the parameter 70 characterizing the broadening of the 
Landau levels in EqlH we have taken 70 = 701/2 be¬ 
cause the broadening of the Eqi transition should have 
contributions from both the n = 0 and n = 1 Landau 
levels. The calculations for the splitting at high mag- 
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netic fields are in excellent agreement with the experi¬ 
mental transmission spectra for both samples S4 (Fig. 
and S5 (Fig. [T]). We neglect /F-phonon absorption [idl) . 
which might account for the discrepancy between theory 
and experiment at the lower frequency and magnetic field 
range of our data. 
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FIG. 9: (color online). Fit of the experimental data (open 
dots) for sample S4 with the conductivity model derived for 
the CDW phase (continuous lines) using Av = 0.24 meV/T 
and Qph = 2.4 meV/T^^^. 


Upon changing the carrier density Ns by a factor of 
~ 10 in our calculations, we find only minor changes 
in the transmission spectra, which is consistent with the 
observation that samples S4 and S5 have remarkably sim¬ 
ilar transmission spectra despite having different carrier 
densities. This is because the disorder induced broaden¬ 
ing reduces the dependence of the sublevel filling factors 
(i^if-f-, etc.) on the carrier density. 

For symmetry breaking driven by electron-electron in¬ 
teractions, the details of the screening function plays 
a vital role in determining the nature of the ground 
state. The additional screening afforded by the multiple 
graphene layers in our epitaxial graphene samples might 
favor the CDW configuration, with two electrons on the 
same sublattice, over the CAF state observed in hBN- 
supported samples Furthermore, coupling between 
rotationally misaligned layers breaks the local A-B sub¬ 


lattice (i.e. valley) symmetry [2^, [s^, promoting the 
CDW ground state. 

VII. CONCLUSIONS 

In conclusion, we have used magneto-optical spec¬ 
troscopy to characterize a SU(4) symmetry broken phase 
in our epitaxial graphene samples. Based on the evo¬ 
lution of the transmission lines near the F phonon fre¬ 
quency, we identify this phase as a CDW phase for the 
specific samples considered, with different occupation 
numbers at valleys K and K'. Because of the valley- 
sensitive nature of the electron-phonon interaction, the 
transmission study used here complements spin-sensitive 
transport measurements in tilted magnetic fields in the 
study of symmetry breaking in graphene. Our experi¬ 
mental method can be applied to open questions such as 
symmetry breaking of the different LLs in graphene and 
bilayer graphene, as well as the effect of disorder on the 
broken symmetry phase in these systems. 
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Appendix 

In general, the detailed analysis of magneto¬ 
transmission spectra requires the use of a multi-layer di¬ 
electric model including all layer dielectric properties of 
the sample. In particular, for each graphene sheet, one 
has to introduce the corresponding components of the op¬ 
tical conductivity tensor tJxx(oj) and axy{uj). Here, the 
X and y-axis lie in the plane of the sample. For instance 
o'xx(<^), in a one-electron approximation, for transitions 
involving the n = 0 LL, is written as: 


^xx .H) 



r,s 


- fsjB)) 

huj — Eps{B) + iTps{B) 


(17) 
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where r,s scan the values 0 and ± 1 , 0 < /,. < 1 is the 
occupation factor of the LL r, M^-^s the optical matrix 
element, = E^ — Eg = Eqi and Trs{B) = 701 (-B) 
measures the broadening of the transition. fx vq 

where is the Fermi velocity given by LDA calcula¬ 
tions In the present work, we have taken for all 

samples vq = 0.85 x 10®ms“^. This is a different param¬ 
eter from vf which appears in Eqi because the energies 
and wavefunctions are corrected to different extents by 
the electron-electron interaction [^. This approach re¬ 
quires the knowledge of the number of effective active 
layers as well as their carrier densities Ng {v = Ng^o/B, 
$0 being the flux quantum) which, in turn, implies some 
approximations. 

The multi-layer dielectric model assumes that each 
graphene sheet is uniformly spread over the sample. This 
is a strong assumption, difficult to justify a priori and 
we have been lead to correct it by assuming a mean cov¬ 
erage which, in the present case for samples S4 and S5, 
has been determined to be about 70 per cent. We next 
evaluate the number Ngff for each sample. In the range 
of magnetic fields 12 to 17 T, the relative transmission 
spectra (Fig. 2, top panel) reaches values above 1 which 
depends on the number Ngf f: we have therefore a guide 
to estimate this quantity. We estimate Neff = 7 for 
samples S4 and S5. 

The carrier density Ng for each layer is determined in 
the following way: one knows that, for 2 < i/ < 6, upon 
increasing B, the intensity of the Eqi absorption starts to 
increase, at the expense of the intensity of the E 12 tran¬ 
sition (£’12 = E 2 — El). The intensity does not change 
with i? for z/ < 2. Therefore, the disappearance of the 
optical transition £12 corresponds to = 2. Following 
the transmission spectra as a function of £, one can eval¬ 
uate the carrier density Wm for each layer m. This is 
an iterative process which converges reasonably (within 
20 per cent) but has to be done independently for each 
sample. The value of Ngi for the layer close to the SiC 
substrate can be set arbitrary to 5 to 6 as given 

by transport data on samples grown under similar con¬ 
ditions: this layer indeed and the two following ones do 
not contribute to the transition £01 in the present ex¬ 
periment. Finally, in the range of magnetic field larger 
than 27 T, where we focus our attention in this paper, 
the number of optically active layers (for optical transi¬ 
tions involving the n = 0 LL) ranges between 3 to 4 for 
samples S4 and S5 with carrier densities ranging from 0.5 
to 12 xlO^^cm-^. 
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